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Abstract 

We present a coarse-grained model for linear polymers with a tunable number of effective atoms (blobs) per 
chain interacting by intra- and inter-molecular potentials obtained at zero density. We show how this model 
is able to accurately reproduce the universal properties of the underlying solution of athermal linear chains 
at various levels of coarse-graining and in a range of chain densities which can be widened by increasing the 
spatial resolution of the multiblob representation, i.e., the number of blobs per chain. The present model is 
unique in its ability to quantitatively predict thermodynamic and large scale structural properties of polymer 
solutions deep in the semidilute regime with a very limited computational effort, overcoming most of the 
problems related to the simulations of semidilute polymer solutions in good solvent conditions. 
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I. INTRODUCTION 



The last decade has witnessed a considerable effort in developing coarse-grained (CG) models 
to bridge the length-scale gap between the microscopic scale and the meso(macro)-scopic scale 
typical of soft-matter and biological systems. Two general strategies have been attempted: the 
structure-based route in which CG models are tuned to reproduce some structural properties at 
the local scale and the thermodynamic-based route in which the models are required to reproduce 
solvation free energies. Also mixed strategies have been developed and applied to a large variety 
of physical systems. An overview of methods and recent applications can be found in Refs. [l-4]. 

The structure-based route groups a number of atoms into effective "particles" and assumes 
state-dependent pair interactions between them. These potentials are derived from the local struc- 
ture of the atomic-level system, using the Iterative Boltzmann Inversion (IBI) , the Inverse 
Monte Carlo (IMC) jfj] or by liquid-state theory techniques 10. 11]. The state dependence of the 



effective pair potentials hides the underlying many-body character and entropic content of the ef- 
fective interactions and poses the question of their transferability. Indeed, if the effective potentials 
are not transferable, setting up the CG model requires deriving the effective interactions for all 
thermodynamic states of interest, seriously limiting the benefits of the coarse-graining strategy. 
Generally speaking, the problem of transferability remains unsolved. 

In the specific case of linear polymers in solutions, structure-based CG models, where a single 
coil is mapped onto a single soft blob, have been introduced quite some time ago 
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18|| . For two 



isolated coils in athermal solvent the effective interaction between the centers of mass of the coils is 
roughly Gaussian, of the order of 2ksT at overlap, and with a range of the order of the coil radius of 
gyration. This single-blob model is only accurate in the dilute regime, in which chain overlaps can 
be neglected. Extension of the single-blob CG model to the semidilute regime can be obtained by 
allowing the use of density dependent pair potentials [3] which can be determined from the radial 
distribution function between the centers of mass of two coils computed from a full monomer 
simulation at finite density. This strategy however suffers for the transferability problem since 
setting up the coarse-grained model requires simulations of the underlying full monomer system 
at all densities of interest. Another limitation of the single-blob model with density dependent 
potential is in representing non-homogeneous systems since the interaction should depend on the 
local density which is not known beforehand and some kind of self-consistent procedure should 
be developed. Furthermore, representing polymers as soft spherically symmetric particles is not 
always appropriate. For instance, in studying polymers adsorbed on surfaces, like polymer brushes 
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or polymer-coated colloids, it is clear that the anchorage to the surface breaks the rotational 
symm etry of the chains, an effect that must be taken into account in any accurate coarse-grained 
model [19j|. A further example is in modelling solutions of A-B block copolymers which cannot be 
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represented as soft particles interacting by a spherically symmetric pair potential 

In a semidilute solution of N linear chains of L monomers in a volume V, chain density, c = N/V, 
is larger than the overlap density c* = 3/4irRg (R g = bL u is the isolated coil radius of gyration, 
b the monomer size and v the scaling exponent) while the monomer density is still very small, a 
condition that can always be satisfied if chains are long enough [2^]. An accurate description of the 
thermodynamic and large-scale structural behavior of polymers in these conditions can in principle 
be obtained by using CG chains, in which a number m of the original monomers are grouped in 
one effective monomer (blob) to map the original chain of L = nm monomers onto a chain of n 
effective blobs (multiblob model). If the level of coarse graining, i.e. the number n of blobs per 
chain, is such that the blob density q, = nc is below the overlap blob density c£ = 3/A7rfg (f g ~ m u 
is the zero-density radius of gyration of the blob), then it is expected that zero-density potentials 
between blobs or different chains can be safely used. These ideas were discussed in Refs. 0- 



301 ] and recently reviewed in Ref. 



311 ] . Although in principle very appealing, the problem with 



this approach is how to obtain the intramolecular potential. Indeed, this inherently many-body 
potential is in principle of increasing complexity when increasing the number of blobs per chain, 
and since blobs are tethered together in some fixed topology, a zero-density expansion cannot be 



311 ] a coarse- 



invoked to decompose it into a sum of two-body, three-body, etc. terms. In Ref. 
grained model for good-solvent conditions with four effective monomers per chain (tetramer) was 
developed and throughly studied. Such a model, set up at zero density, was found to be accurate 
up to reduced polymer density $ = c/c* ~ 2 supporting the multiblob ideas. 

In this paper we present the extension of the tetramer model to chains with an arbitrary number 
of blobs, an extension needed to explore the semidilute regime at high polymer concentration. We 
show that a careful parametrization of the n-body interactions up to n = 4 is enough to obtain a 
very accurate CG potential field fully transferable both with the number of blobs per chain and 
with the density of chains, provided that the length of the effective chains is chosen to ensure that 
the reduced blob density remains always low. The paper is organized as follows. In section II we 
will provide a theoretical framework on which any accurate coarse-grained model should be based. 
In the following section III we present our multiblob model and report results for systems of chains 
of varying length both at zero density (section III. a) and at finite density in the semidilute regime 
(section Ill.b). In section IV we collect our conclusion and perspectives. Finally an appendix 
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reports the tetramer potentials as obtained by the IBI and their explicit parametrization. 



II. THEORETICAL FRAMEWORK 

A crucial requirement of any accurate coarse-grained representation (CGR) representation of 
the underlying full-monomer (FM) chain model is to preserve the value of the radius of gyration 
for all values of the number of blobs n. Let us consider a chain of L monomers mapped onto 
a chain of n blobs, each representing the center of mass of a subchain of m monomers. The 
fundamental relation R 2 g = R^ b (n) + r^(n) among the chain radius of gyration R g , the radius of 
gyration R g ^{n) of the chain of n blobs, and the mean blob radius of gyration r g (n), holds for any 
single configuration and therefore on average. Here rz{n) is the average of the square radius of 
gyration of all blobs in the chain 

1 n 
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where rj is the position of the i-th monomer in the chain. If Sj is the center of mass of the i-th 
blob in the chain, 
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Rg,b(n) is defined by 



In the scaling limit, it has been found that [3l| 
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where ^ = 0.587597(7) 32|. Here and in the following we will use a hat to indicate zero-density 



averages. For n » 1, R 9t b(n) ~ fi 9 and is! 2 6 (n)/f 2 (n) = 0.94n 2,y . 

Any consistent multiblob (MB) model must obey these relations when varying n. A proper def- 
inition of R g , independent on n, including the prefactor, is crucial since chain structural properties 



are universal only if distances are expressed in terms of rescaled distances p = r/R g . Moreover, 
at finite density results from different values of n should be compared at the same value of the 
polymer volume fraction $ = c/c* = (4TiRg)c/3. Changing the definition of R g with n, changes the 
definition of both p and <3? and the comparison among models with different n becomes meaningless. 

In Ref. 3l|] a MB model with four blobs per chain (tetramer) was developed and throughly 
studied. Central potentials between first, second and third neighbors along the chain, as well as a 
bending and torsional angle potentials were determined by IBI to reproduce the pair distances and 



ail. 



angular distributions of an isolated FM chain in the scaling limit, mapped onto four blobs 
Moreover, an intermolecular Gaussian potential between any pair of blobs of different tetramers 
was assumed and optimized to reproduce the radial distribution function between the centers of 
mass of two isolated chains in the scaling limit 17|] . The model perfectly reproduces the structure 



of the chains and provides accurate results for the thermodynamics of the solution up to <3? ~ 2. 
The non-uniform (universal) angular distributions of the tetramer are induced by the specific 
mapping of subchains onto their centers of mass, a procedure which introduces an explicit angular 



correlation along the chain even for ideal chains 



Failure in reproducing the correct angular 



distributions is at the basis of the inconsistent behavior observed for simpler multiblob models 



28, 



29|, 



311 ] . Note that while central potentials are two-body interactions, bending and torsional 



ang le potentials represent genuine three- and four-body interactions, respectively. Interestingly 



311 ]. a hierarchical order in the intensity of the various potentials is observed, the strongest and 
most relevant being the central first-neighbor interaction (bonding), followed by the second- and 
third-neighbor central interactions, and finally by the bending and the torsional angle potentials. 
Therefore, correlations along the chain decrease with the chemical distance and many-body effects 
are smaller than two-body ones, although of different nature. This observation suggests to use 
tetramers as building blocks of a resolution-invariant multiblob model. 

III. TETRAMER-BASED MULTIBLOB MODELS AND RESULTS 
A. Zero-density transferability 

A first tetramer-based multiblob model (4MB-1) is obtained by using the tetramer potentials 



of Ref. 3l|] reported in the Appendix. In model 4MB-1 one considers six different intramolecular 
potentials, each of them depending on a single scalar variable. We use four bonding pair potentials: 
blobs i and j interact with a central potential Vij(p;4) with p = |s$ — Sj\ /R g . Symmetry imposes 
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^i,2(p;4) = V3,4(p;4) and V\^{p;4) = V2,4(/o;4) so that there are only four different central poten- 
tials. Then we use a three-body potential Vf,(cos/3j) acting on the two equivalent bending angles 
and a four-body potential Vt(0) acting on the torsional angle, where /3 and 9 are defined as 

cos & = "ifT^ (2 = 1,2), (7) 

cqs9 = (bix^.Q^xba) 
|bi x b 2 ||b 2 x b 3 | 

where bj = Sj+i — Sj is the bond vector between blobs i ad i + 1. The intermolecular potential 
is instead represented by a unique blob-blob pair interaction W(p; 4) acting between any pair of 
blobs of different tetramers. 

When using this potential field to build chains of n blobs, the blob radius of gyration r g is taken 
as the length scale of the central potentials. In practice this means that when transferring the 
potentials from the model with 4 blobs to the model with n blobs, the following identification is 
adopted 

V ij (p;n) = V ij (p n ;4), (9) 
W(p;n) = W( Pn ;4) (10) 

where, according to Eq. (JB]), p n is defined by 

/4V 1.03-0.04/4 
Pn=P {n) 1.03 - 0.04/n (U) 

The angular potentials are instead invariant under this scale transformation. 

This model, 4MB-1, is not found to be fully satisfactory, as illustrated in tablelH where we report 
results for the zero-density chain radius of gyration and the second dimensionless virial coefficient 



for increasing n for polymer chains in the scaling limit (FM) [341 ] . the multiblob model (4MB-1), 
and the modified multiblob model (4MB-2) still to be discussed. In particular, in model 4MB-1, 
Rg.b( n ) I Rg becomes greater than 1 for increasing n, a fact clearly incompatible with the multiblob 
procedure [see Eq. (J5])]. Also a systematic increase with n in the second virial coefficient of 4MB-1 
is observed which witnesses the inability of such model to be fully consistent at all resolutions. 

The origin of such inconsistency can be ascribed to a residual inaccuracy of the original tetramer 
model to reproduce four-body correlations. In Fig. [1] we show, for n = 4, the probability distri- 
bution of cos /?i3 = (bi • t>3) /(|t»i | |t»3 1). Data of the 4MB-1 model are compared with FM-CGR 
results. Despite the explicit presence of the angular potential which reproduces very accurately 
the torsional angle distribution 3l|], P(cos/3i3) is not well reproduced in the 4MB-1 model — more 
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TABLE I: R g ^(n)/Rg and universal ratio A 2 = B 2 jB? g (B 2 is the second virial coefficient) as a function of 
n for polymers obtained by full-monomer simulations (FM-CGR), and for the original and modified 4MB 
models. The universal value for polymers of the second virial coefficient combination is A 2 = 5.500(3) 



A 2 = B 2 B 3 „ 



n FM-CGR 4MB-1 4MB-2 ^ 2 (4MB-1) A 2 (4MB-2) 
4 0.8921(1) 0.8933(3) 0.8903(2) 5.597(1) 5.596(3) 
10 0.9642(1) 0.9812(3) 0.9638(7) 5.676(2) 5.619(2) 
20 0.9842(1) 1.0060(6) 0.983(3) 5.714(2) 5.614(3) 
30 0.9902(1) 1.0131(8) 0.990(3) 5.734(3) 5.618(6) 
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FIG. 1: Probability distribution of cos/?i3 in the original tetramcr model (4MB-1) and in the modified 
tetramer model (4MB-2). The distribution for a CGR of a full-monomer chain with four blobs is also 
reported for comparison (FM-CGR). The inset displays the relative difference of the two models with the 
FM-CGR prediction. 



elongated configurations are apparently enhanced. This inaccuracy, irrelevant for the properties 
of the multiblob model at the tetramer level (re = 4), accumulates when transferring the tetramer 
potentials to longer chains, producing the observed increase of R gt b(n)/R g and A2 with n. A sec- 
ond, modified tetramer model which reproduces the FM-CGR behavior can be obtained by adding 
a new potential on cos/?i3. As for the other potentials, the IBI procedure is applied to extract 
the new optimal potential. At the same time the dihedral angle potential is re-optimized to keep 
the level of accuracy of the original tetramer model. Explicit expressions for the potentials are 
provided in the Appendix. The good accuracy of the new model is illustrated in figure [TJ The 
modified tetramer model can now be safely used as the building block of a fully consistent model 
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(4MB-2) which puts in action the ideas of the multiblob approach. As illustrated in table HI at 
zero density the agreement with the FM-CGR predictions for the rescaled radius of gyration is 
excellent for all n values investigated. Also the virial coefficient is independent on n and in good 
agreement with the FM prediction A2 = 5.500(3) 36]. This is the first requirement of transferabil- 



ity with n at zero density. Note that the third virial coefficient appears to be much less sensible to 
model inaccuracies: we obtain A3 = RgBj = 10.3(2) both for 4MB-1 and 4MB-2 for all n values 



36J: A 3 = 9.80(2). 



investigated to be compared with the FM prediction 

To check that our parametrization exhausts all four-body terms we compare in Fig. [2] the 
cross correlation between the cosine of the two bending angles of the tetramer as obtained by the 
FM-CGR and by the 4MB-2 model. This quantity is defined as 

(5(cos f3 ± - x)5(cos p 2 ~ y)} (m 

In the absence of correlation we have C(x, y) = 1. We observe the presence of a region of depleted 
correlation near f3\ = P2 = 0, and two regions of enhanced correlation around f3\ = 0, $2 = 7T and 
Pi = tx^2 = 0, respectively. The 4MB-2 model reproduces fairly well the results of the FM-CGR 
(see the right panel), except in a very narrow region around f3\ = = where statistical noise is 
also larger, indicating the accuracy of our parametrization. 




FIG. 2: Cross-correlation function C n (x,y). Comparison between the 4MB-2 model (left panel) FM-CGR 
chains (central panel) for n = 4 at $ = 0. In the right panel, we show the relative deviations of the 4MB-2 
results from the FM-CGR ones, AC{x,y) = 100\C 4 MB-2(x,y)/C FM (x,y) - 1|. 



To further investigate if of the new model 4MB-2 is able to reproduce FM-CGR correlations, 
we computed genuine five-body intramolecular correlations for varying number of blobs n. Note 
that two-, three- and four-body correlations at short chemical distance are explicitly built into the 
model, and five-body correlations are the first missing terms. A useful tool to characterize these 
correlations is the torsional angle correlation function % n ($i, $2) between two subsequent dihedral 



angles along the backbone. If B{ is the i-th torsional angle along the chain, we adopt the following 
definition for the function H n (&i, $2)'- 



S 4 (^l-^))W*2-%)> 



where (...) is the average over single-chain configurations, and the subscript n emphasizes the 
dependence of T~i n on the number n of blobs. It is possible to define similar functions for any 
pair of torsional angles along the backbone, not necessarily nearest neighbors. Since we expect 
the correlation between such pairs to decay rather quickly with the chemical distance we limit our 
analysis to nearest neighbors. Note that with this definition, 7i n (^i, $2) = 1 in the absence of 
correlation. In Fig. [3] we report the function T~L n for n = 10, 20, 30 in both 4MB-2 (left panels) and 
FM (central panels) representation and the relative deviations (right panels). The 4MB-2 model 
predictions of H n are quite close to the polymer results, the deviations observed are at most 6% in 
well localized regions of the (<J>i, $2) plane, and are nearly independent on n, confirming the good 
transferability of the potential, as observed for the other structural properties. This observation 
supports our implicit assumption that irreducible n-body terms in the intramolecular potential for 
n > 5 are negligible. 

It might appear surprising that we spent so much effort in modeling the intramolecular potential 
while the intermolecular interaction is treated in a very simple way. The effective interaction among 
N chains, of n blobs each, is a general function of (3nN — 6) variables. According to the multiblob 
ideas, at any finite reduced chain density <3>, we can adopt the small density expansion and neglect 
many-chain interaction terms for n large enough. However the effective intermolecular potential 
between two chains of n blobs remains a function of 6n — 6 (n > 1) variables. In our model 
we adopted several simplifying hypotheses: we replace the intermolecular chain-chain interaction 
by a sum of pair- wise potentials between blobs of different chains, we assume that all these pair 
potentials can be adequately approximated by a single central potential, and finally we approximate 
the potential with a single Gaussian function. It would not be too difficult to avoid the last two 
assumptions. It is instead crucial, from a practical point of view, to assume pairwise blob-blob 
interactions. This unique intermolecular blob-blob potential is optimized in such a way to reproduce 
the polymer center-of-mass radial distribution function (RDF) gcM determined from full monomer 



simulations 



311 ] . This criterion ensures the correct thermodynamics for the CG model since the 



isothermal compressibility is related to the integral of the RDF by the compressibility equation 



3J 



(^y 1 = 1 + £;*J(90M(p)-i)*p. (i4) 
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FIG. 3: W n ($i,$ 2 ). Comparison between the 4MB-2 model (left column) and FM-CGR chains (central 
column) for n = 10 (upper row), n = 20 (central row) and n — 30 (lower row) at $ = 0. On the 
right column, for each value of n, we show the deviations of the 4MB-2 results from the polymer ones, 
A($i,$ 2 ) = 100|H 4 mb-2($i,$2)/Hfm($i,$2) - 1| 



where IT is the osmotic pressure of the solution and /3 = 1/ksT. An equally good choice as 
thermodynamics is concerned, would have been to reproduce the RDF between any pair of blobs 
of different chains since in Eq. (|14p . gcM can be replaced by any other intermolecular RDF 381 ] . 
With our simplifying hypothesis of a unique blob-blob intermolecular interaction we are of course 
unable to perfectly reproduce RDFs between center-of-mass and between any pair of blobs of the 
two chains simultaneously. In figure 0] we show some of those RDFs predicted by the 4MB-2 model 
and we compare them with the corresponding quantities from FM simulation. We observe perfect 
agreement for the RDF between the centers of mass but, as expected, some small deviations 
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occur for the other RDFs at distances p < 0.5, a range irrelevant to thermodynamics. These 
small deviations are the price to pay in order to have a much simplified model, still accurate as 
thermodynamics is concerned. 
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FIG. 4: Radial distribution functions (RDF) between a pair of isolated chains (zero density) as a function of 
p = r/Rg. Comparison between FM-CGR and 4MB-2 model predictions for n = 4 (left panel) and n = 10 
(right panel). RDF between centers of mass (CM-CM), between end monomers (E-E), between central 
monomers (C-C) and between end and central monomers (E-C) are shown. 



B. Model transferability with density 

The second and more difficult requirement for a transferable MB model is the ability to re- 
produce FM predictions at any finite reduced chain density by increasing the number n of blobs. 
To illustrate the ability of our model to match this second requirement we compare single-chain 
properties and the Equation of State (EOS) with FM predictions in the reduced density range 
$ = c/c* < 9. Note that for the present MB model, which is resolution invariant at zero density, 
the definition of the reduced polymer density <I> is given in terms of an n-independent characteristic 
length, the FM radius of gyration R g . This implies that at given values of MB models with 
different n are all at the same density c = N/V of the underlying polymer model. Conversely, 



311 ]. comparing systems with different n at the 



when adopting non-resolution invariant models [28, 
same value of implies comparing systems at different values of the absolute polymer density c 
which is inconsistent. 

In Fig. [5]we compare FM-CGR and 4MB-2 predictions for the intramolecular radial distribution 
function p 2 gi n tra{p) , which is universal in terms of the reduced intra-blob distance p = r/R g , for 
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FIG. 5: Intramolecular radial distribution function p 2 gi n tra{p) versus the reduced blob distance p — r/R g 
for model 4MB-2 with n = 20 and for the FM-CGR chain with 20 blobs. Four values of the reduced density 
are reported: $ = 0,1.09,4.36,8.72. For sake of clarity, results at different densities are shifted upward 
according to the legend. 



a chain of n = 20 blobs and for increasing density. The observed agreement is remarkable in 
particular at the highest density investigated, $ = 8.72, deep inside the semidilute regime. The 
same kind of agreement is observed for the distribution of the radius of gyration P(R/R g ) shown in 
Fig. [6l The agreement worsens at the highest densities for shorter chains, but it remains very good 
for n > 20. We would like to emphasize the nontrivial nature of the observed agreement. Indeed, 
P 2 Qintra{p) is a weighted average of the length distribution over all pairs of blobs, so the agreement 
for this quantity demonstrates the accuracy of the 4MB-2 model to reproduce pair distances well 
beyond four neighbors. P(R/R g ) includes information on pair distances but also on all possible 
cross-terms between different pairs, so the observed agreement for this quantity is even less obvious. 
Finally, we have computed the EOS, using the standard molecular virial route to the pressure 



for several values of n 



391 ]. In Table HI] we report the compressibility factor Z = j3H/c, where II 



is the osmotic pressure of the solution. Up to density <I> ~ 2, results for all values of n studied 
differ by less than 1% from the asymptotic polymer value. When density increases longer chains 
are systematically more accurate. At $ = 4.36 the relative deviation in Z is ~ 5% for n = 10, but 
< 1% for both n = 20 and n = 30. At $ = 8.72, the deviation is ~ 15% for n = 10, and decreases 
to ~ 10% for n = 20 and ~ 6% for n = 30. This is expected from the multiblob heuristic argument 
since from Eq. ([6]) the reduced blob density, 
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FIG. 6: Chain radius-of-gyration probability distribution function P(R/R g ) (R is the radius of gyration 
of a single chain and R g its zero-density average) at $ = 0,1.09,4.36,8.72 for a chain of n = 20 blobs. 
Comparison between the 4MB-2 model and the FM-CGR predictions. 



TABLE II: Compressibility factor Z for the 4MB-2 model for several values of n and for polymers in the 
scaling limit 



at various reduced densities $ from the dilute to the semidilutc regime. 



n = 10 n = 20 n = 30 



FM 



0.54 1.875(1) 1.874(2) 1.870(2) 1.85(1) 

1.09 2.9760(4) 2.9784(6) 2.9819(1) 2.96(1) 

2.18 5.6004(7) 5.708(1) 5.729(3) 5.63(1) 

4.36 11.4541(5) 12.058(1) 12.269(3) 12.3(1) 

6.54 17.532(2) 18.931(4) 19.515(5) 20.0(1) 

8.72 23.751(1) 25.592(4) 26.989(8) 28.7(1) 



decreases, at fixed for increasing n as long as v > 1/3. At $ = 8.72, ^(lO) ~ 1.6, ^(20) ~ 1.0 
and <J>b(30) ~ 0.7 which explains why increasingly large n-values are necessary in order to neglect 
three- (or more)-blobs contacts as in the present model. 

It is of some interest to investigate the relation between the blob radius of gyration of the 
coarse grained representation r g (n) and the density-dependent correlation length £ characteristic 



of semidilute polymer solutions 



271 ] . also called correlation blob size. The latter quantity can be 



inferred from the single chain structure factor, 5^ n t ^0 (q , ) = L e jq riJ ), as the characteristic 

length at which a crossover from the ideal behavior (~ q 2 ) to the excluded- volume behavior 
(~ q~ l / v ) is observed j^oj-Q- I n figure [7] we show at <I> = 8.72 the Kratky plot, (qR g ) 2 Si n tra(q) 
against qR g (&), for the full monomer Domb- Joyce model (FM), the FM-CGR with n = 30 and 
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FIG. 7: (color online) Kratky plot at <f> = 8.72 as a function of Q = qR g (&). Full monomer Domb- 
Joyce model (FM, red squares), the same chain model mapped onto a 30-blobs chain (FM-CGR, blue line), 
the model 4MB-2 with n = 30 blobs (open green circles), the Debye function (dotted line). Note that 
Q = qR g .b{&) has been used for the FM-CGR and 4MB-2 models. For n = 30 the difference between 
between R g and R g ,b is already small and will become even smaller for larger values of n. 

the 4MB-2 model with the same number of blobs. In the FM-CGR case Si n t r a{o) is defined as 
above with monomer positions replaced by the blob center-of-mass positions Sj. The correlation 
blob size £ can be readily determined from the FM data to be £/fL(8.72) ~ 2-7r/10 ~ 0.63, which 
corresponds to £,/R g — 0.53 since R g (8.72) / R g = 0.8423 at this reduced density. On the other hand 
from Eq. ([6]) we obtain f g {n)/R g = k/n u = 0.27, 0.18, and 0.14 for n = 10, 20, and 30, respectively, 
which are considerably smaller than £. This result is consistent with the multiblob argument: 
:'rom Eq. (I15p . <&b oc n^ 1-3 ^^ = (R g /r g )^ 3u ^ u ^. On the other hand in the semidilute regime 
3] i/R g oc <3W( 1_3l/ ). Substituting the former relation in the latter we get: ® b oc {r g /C)^ u ~ l ^ h '■ 
Therefore, as far as v > 1/3, 1 implies r g <C £ in agreement with our findings. 

In Fig. [7] we note that the CGR at such resolution reproduces FM data up to Q ~ 2. If one 
wishes to have access to more local properties he should take a CGR with more blobs. In order 
to follow the crossover from ideal to excluded volume behavior in the CGR a quite larger number 
of blobs should be considered. Note that the large-g behavior of the CGR (~ q 2 ) is universal 
since our blobs contain a sufficiently large number of monomers to be representative of the scaling 
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limit. Therefore this large-g behavior is inherent to the CGR and not to the particular MB model 
we considered, namely the 4MB-2 model. Finally, note the that perfect agreement between the 
4MB-2 prediction and the FM-CGR data for the single chain structure factor is fully compatible 
with previous results on gi n tra(r) (see Fig. [5]). 



IV. CONCLUSIONS 



In conclusion, we have developed a fully consistent, scale preserving, multiblob model for linear 
polymers in good solvent. The model is built by transferring to smaller length scales a tetramer 
model parametrized to reproduce a number of scalar correlations of the FM chain at zero density. 
We have shown that this model is fully consistent when varying the number of blobs at zero 
density, and more relevant, it is able to reproduce the universal EOS for athermal semidilute 
polymer solutions at high chain concentrations if a sufficiently large numbers of blobs per chain 
is chosen in such a way to always remain in the dilute blob regime. In particular, the present 
multiblob model with only 30 blobs provides the compressibility factor at ~ 9 with an accuracy 
of ~ 5%, a level of accuracy which would require the use of the order of thousand monomers even 



with the most efficient lattice model, name 
the leading-order corrections to scaling 



y the Domb-Joyce model tuned in such a way to cancel 
351 ] . The present CG strategy may be extended to more 
complex systems and thus opens the way to a quantitative study of semidilute polymer solutions in 
situations where local structure is important like, for instance, colloid-polymer mixtures, polymer- 
decorated colloidal systems and polymer brushes at relatively low grafting density. It would also 
be extremely interesting to develop an analogous multiblob strategy to describe the cross-over 
between the good solvent and the 9 solvent conditions and finally to extend this approach to 
diblock copolymer solutions. Work in these directions is in progress. 

As a final remark it is interesting to observe that our model represents, to a very good approx- 
imation, the first example of a truly "fixed-point" model in the renormalization-group language, 
the elusive model which reproduces the scaling behavior at any level of coarse-graining. Previous 
work focused on defining models in which the leading scaling corrections were zero within errors. 
Here, instead, we obtain a model in which all scaling corrections are apparently small, and which 
is thus able to predict the scaling behavior for any n > 4. Of course, there is a price to pay: for 
any n, we do not have access to all possible observables, but only to the large-scale properties that 
can be modelled by the chosen CG representation. 
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VI. APPENDIX: EXPLICIT EXPRESSIONS OF THE TETRAMER POTENTIALS 

In this appendix we provide full details on the potentials of the two models used in the present 
work, namely 4MB-1 and 4MB-2. Distance in the central potentials is p = r/R g . 

The bonding potentials /3Vi2(p; 4) = fiVz^p; 4) and pV2z{p; 4) have been parametrized for p < 3 

as 

PVij(p; 4) = (V -a - m) exp(-a 2y o 2 ) + ai exp(-a 3 p 2 ) + a 4/ o 2 + a 5/ o 4 + a , (16) 

where Vo is the value of the potential at the origin, ao is a constant which has been fixed so that 
the potential vanishes at the minimum, and we only consider the range p < 3. 

TABLE III: Numerical coefficients for the tetramer potentials (n = 4) for 4MB-1 and 4MB-2 models 



Model 


Potential 


V 


a 


ai 


"2 


03 


a 4 


as 


4MB-l,4MB-2 


v 23 


0.8808 


-6.8935 


1.8693 


0.4287 


5.3453 


4.4274 


-0.0202 


4MB-l,4MB-2 


Vl2 


1.0392 


-3.9339 


1.54603 


0.7602 


6.8442 


3.3646 


-0.0194 


4MB-l,4MB-2 


Viz 






-3.9518 


2.9224 


-1.1555 


0.1321 




4MB-l,4MB-2 


U 






0.2238 


-0.2284 


0.0769 


0.0600 




4MB-1 


Vi 




0.1425 


1.0698 


4.5475 


0.0391 


0.14345 




4MB-2 


Vi 




-0.0797 


0.970 


0.990 


-0.0251 


0.0785 




4MB-2 


v bl3 






0.1593 


-0.0371 




-0.0158 





The potentials Vis(p;A) = V24(/o;4) between next-to- nearest neighbors has been parametrized 
as (p < 3) 

pVi 3 (fi\ 4) = PV 2 i{p; 4) = aip + a 2 p 2 + a 3 p 3 + a 4 p 4 . (17) 
The bending potential is parametrized as 

/3I4(cos (3) = ai(cos j3 — a 2 ) 2 + a3(cos /3 — a 2 ) 6 + a4(cos j3 — a 2 ) 7 (18) 



16 




FIG. 8: Potentials for the tetramer 4MB-2 and 4MB-1. In the left column we show the pair potentials 
between first-neighbors along the chain (top) and between ncxt-to- nearest neighbors (bottom) , as a function 
of p = r/Rg. In the central column we report the bending-angle potential, the potential acting on the 
cos^Ais) as a function of cos /? (top), and the torsion-angle potential as a function of $ (in radians) for both 
the models 4MB- 1,4MB- 2 (bottom). In the right column we report the intermolecular pair potential (top) 
and the potential between the first and the last atom of the tetramer (bottom) as a function of p = r/R g . The 
symbols represent the numerical results obtained through IBI procedure and the lines are the corresponding 
interpolations. 

(the minimum of the potential corresponds to cos/3 = 0.2), while the torsion potential is given by 

PV t {®) = a (l + 03$) sin(ai$ + a 2 ) + a 4 . (19) 

The potential acting on the cos/3i3 has been parametrized as: 

PV bl3 (cos /3 13 ) = m cos /3 i3 + a 2 cos 2 /3 i3 + a 4 . (20) 

Numerical values for the coefficients for n = 4 are reported in Table IIIII and the potential are 
displayed in figure [8l 
Finally, 

pVu(p;4) = 1.86exp(-4.08425p 2 ), /3W(p; 4) = 1.66 exp(-3.9p 2 ). (21) 
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